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SUMMARY 

A  low  speed  investigation  of  the  flow  over  aspect  ratio  one  delta  wings  of 
varying  thickness  has  been  made  to  better  understand  the  relation  between  the 
vortices  produced  by  leading  edge  separation  and  the  non-linearity  of  the  lift 
curve.   The  formation  and  the  position  of  the  vortex  cores  were  determined  by 
means  of  smoke  flow  visualization  techniques.   Lift  curves  of  the  models  were 
obtained  from  the  wind  tunnel. 

It  was  found  that  at  the  same  angle  of  attack  the  vortex  cores  moved  outward 
on  the  wing  as  wing  thickness  increased.   Downstream  of  the  trailing  edge,  the 
vortex  cores  followed  a  helical  path  which  is  believed  to  be  related  to  the 
"vortex  explosion"  phenomenon.   An  empirical  equation  was  developed  which  predicts 
the  lift  curves  for  sharp  leading  edged  delta  wings  of  various  thicknesses  and 
of  aspect  ratios  from  one  to  two. 


INTRODUCTION 

Traditionally,  the  influence  of  viscosity  on  the  flow  over  a  wing  was  assumed 
to  be  confined  to  a  thin  boundary  layer.   Outside  the  boundary  layer,  the  flow  was 
considered  to  be  inviscid  and  could  be  analyzed  by  linear  potential  theory. 
Where  the  leading  edge  is  sharp  on  a  thin  swept  wing,  the  flow  may  not  be 
considered  inviscid,  because  the  boundary  layers  developed  on  the  lower  surface 
separate  at  the  sharp  leading  edge.   Specifically,  for  a  delta  wing,  with  highly 
swept,  sharp  leading  edges,  at  angles  of  attack  greater  than  zero,  the  boundary 
layers  on  the  pressure  side  separate  at  the  sharp  leading  edges  to  form  vortices 
which  grow  in  size  and  strength  toward  the  trailing  edge.   The  leading  edge 
vortices  are  the  significant  characteristic  of  such  wings. 

The  effect  of  the  vortices  is  to  produce  a  complex  flow  field  over  the  wing 
with  associated  non-linear  aerodynamic  characteristics.   For  such  a  flow  field, 
deviation  from  potential  wing  theory  exists  which  results  in  a  non-linear  lift  curve, 

Following  the  principle  that  only  when  the  physical  situation  is  fully 
understood  can  the  mathematical  treatment  be  applied  intelligently,  smoke 
visualization  techniques  were  applied  to  the  flow  field.   In  addition  to  the 
qualitative  photographs,  quantitative  results  were  obtained  in  the  form  of  vortex 
core  locations  with  angles  of  attack  for  the  three  aspect  ratio  one  delta  wings 
shown  in  F  ig.  I . 

Based  on  the  vortex  core  locations,  an  empirical  expression  for  the  lift 
curve  was  obtained  which  is  applicable  to  delta  wings  up  to  12  percent  thick  and 
of  aspect  ratio  one  to  two. 


2. 

< 

BACKGROUND 
H  i story 

In  1936,  Winter   first  described  the  conical  vortex  at  the  wing  leading  edge. 
By  means  of  flow  visualization  techniques  (smoke,  tufts,  and  "soot  and  petroleum"), 
he  described  the  flow  as  "flowing  spirally,  forming  a 'vortex  braid'  whose  position, 
magnitude  and  extent  are  determined  by  the  plan  form,  aspect  ratio  and  angle  of 

attack." 

2 
In  1947,  Wilson  and  Lovell   determined  by  tests  of  a  delta  wing  and  tuft 

analysis,  that  "the  flow  over  triangular  wings  of  low  aspect  ratio  is 

characterized  by  vortices  above  the  upper  surface  of  the  wing,  inboard  of  the  tips. 

These  vortices  aid  in  obtaining  high* maximum  lift  coefficients  and  can  be  produced 

by  using  airfoil  sections  having  sharp  leading  edges." 

3 
Ornberg   in  1954,  identified  the  secondary  vortices  produced  by  the 

separation  of  the  upper  surface  boundary  layer. 

Many  have  contributed  significantly  to  the  understanding  of  the  flow,  in  the 

form  of  wind  tunnel  tests,  flow  visualization  techniques  and  mathematical  treatment, 

A  partial  listing  of  the  papers  published  on  the  subject  is  contained  in  the 

references  and  bibliography  section. 

i". ■■..'•-.  .       : 
Flow  Description 

As  the  angle  of  attack  increases  from  zero,  the  stagnation  point  moves  from 

the  leading  edge  to  the  under  surface  of  the  wing.   The  flow  proceeds  from  this 

stagnation  point  to  the  leading  edge  due  to  the  favorable  pressure  gradient  in 

this  region.   Suction  on  the  upper  surface  tends  to  draw  the  flow  around  the 

leading  edge.   However,  for  a  sharp  leading  edge,  the  centrifugal  force  on  a  fluid 

element  approaches  infinity  and  the  fluid  cannot  remain  attached  to  the  surface 

of  the  wing.   The  separated  fluid  takes  the  form  of  a  vortex  sheet  which  grows  in 

strength  as  it  proceeds  downstream  on  the  body.   Thus,  when  it  rolls  up  it  takes 


the  form  of  a  conical  spiral.   The  vortex  sheets  roll  up  in  opposite  senses  on 
the  wing  into  cores  of  rotating  fluid  which  lie  above  and  inboard  of  the  leading 
edges  as  seen  in  Fig.  7.   These  leading  edge  vortices  are  the  dominant  feature 
of  the  flow. 

In  addition  to  the  strong  flow  from  the  lower  surface,  there  is  the  secondary 
flow  pattern  which  is  produced  by  the  separation  of  the  outflowing  air  in  the 
boundary  layer  on  the  upper  surface  as  it  nears  the  leading  edge,  as  seen  in 
Fig.  2.   This  separation  occurs  due  to  an  adverse  pressure  gradient  on  the  upper 
surface  outboard  of  the  primary  vortex  cores.   The  secondary  separations  give  rise 

to  small  triangular-shaped  viscous  regions  on  the  upper  surface  near  the  leading 

4 
edge.   Investigation  of  these  viscous  regions  by  means  of  a  pressure  probe 

located  the  secondary  vortex  cores  below  and  outboard  of  the  primary  cores  and 

confirmed  the  rotation  of  the  vortices  as  opposite  in  sense  to  the  primary 

vortices.   Downstream  of  the  trailing  edge,  the  secondary  vortices  are  well 

defined,  as  seen  in  Fig.  II. 

Theoretical  Treatment 

The  flow  pattern  in  quite  different  from  that  usually  considered  in  the 

5 
theory  of  inviscid  fluid  motion  presented  by  R.  T.  Jones.    Jones'  theory  requires 

the  fluid  to  turn  through  an  angle  at  a  sharp  leading  edge  and  so,  in  the 

linearized  approximation,  predicts  infinite  values  of  the  velocity  and  pressure 

along  it.   It  is  not  these  singularities  in  themselves  which  render  such  theories 

for  the  slender  delta  wing  inadequate,  since  they  are  of  the  type  acceptable 

elsewhere  in  the  theory  of  thin  airfoils;  but  rather  that  the  flow  pattern  found 

in  a  real  fluid  does  not  resemble  that  of  the  potential  flow  model  which  predicts 

them. 


4. 

Mangier  and  Smith  treat  the  flow  theoretically  by  constructing  a  potential 
flow  model,  in  which  the  vortex  layer  is  replaced  by  a  vortex  sheet  of  spiral 
form  and  the  problem  is  then  reduced  to  a  two—  d  imens iona I  one  by  the  use  of 
slender  body  theory  and  the  assumption  of  a  conical  vorticity  field.   The  shape 
and  the  strength  of  the  vortex  sheet  are  determined  by  utilizing  two  boundary 
conditions:  the  sheet  lies  in  a  stream  surface  of  the  flow  with  no  difference 
in  the  pressure  across  it  and  the  velocity  is  assumed  to  be  constant  along 
strraiight  lines  through  the  apex  of  the  delta.   This  assumption  agrees  with 
experiments  at  supersonic  speeds,  however,  at  low  speeds,  this  assumption  breaks 
down  near  the  trailing  edge  as  indicated  in  Fig.  33.   It  is  noted  that  recent 
tests  at  the  National  Physical  Laboratory  (Aero.  Div.)  have  found  core  velocities 
in  the  direction  of  the  core  path  to  be  on  the  order  of  1.5  times  the  free  stream 
velocity. 

o 

Brown  and  Michael   considered  the  actual  flow  solution  to  be  extremely 
difficult,  hence  a  flow  model  more  amenable  to  calculation  was  adopted.   This 
simplified  model  replaced  the  spiral  sheet  with  two  concentrated  line  vortices 
above  the  wing  and  two  feeding  vortex  sheets  connecting  the  source  of  vorticity 
(leading  edge)  and  the  concentrated  line  vortices.   Here  also,  conical  flow  is 

assumed. 

9 
Legendre   in  1952  considered  the  problem  as  a  potential  flow  and  represented 

the  vorticity  in  the  fluid  by  a  pair  of  isolated  vortices  lying  along  streamlines. 

This  introduced  a  many-valued  function  for  the  pressure  which  placed  a  pressure 

discontinuity  along  the  path  around  either  one  of  the  isolated  vortices.   In 

10 
1953,  Adams   suggested  that  these  discontinuities  should  be  placed  on  surfaces 

or  sheets  joining  the  vortices  to  the  leading  edge.   Following  Adams'  suggestion, 

Legendre  published  a  revised  analysis  of  the  problem  in  1953. 


5. 

Figure  3  shows  the  predicted  inward  and  upward  movement  of  the  primary  vortex 
cores  with  angle  of  attack  for  slender  delta  wings  with  leading  edge  separation. 

The  theoretical  curves  are  derived  for  supersonic  flow  with  the  secondary  vortices 

I  I 
not  considered.   Experimental  curves  from  Fink  and  Taylor   and  this  report  are 

included  for  comparison.   It  is  noted  that  the  experimental  curves  are  from  low 

speed  studies. 

Authoritative  opinion  on  the  effect  of  the  secondary  vortices  is  quite 

diverse.   Reference  4  states  that  the  real  flow  was  different  from  that  derived 

1 2 

by  Mangier  and  Smith, due  to  secondary  vortices.   Lee   suggests  that  the  influence 

on  the  main  flow  by  secondary  vortices  is  negligible.   Attention  is  directed  to 
Figs.  9  and  II  which  show  the  secondary  vortices  downstream  to  be  of  comparable 
size  to  the  primary  vortices.   This  indicates  that  the  size  and  hence  the  strength 
of  the  secondary  vortices  on  the  body  may  be  of  appreciable  magnitude  and  hence 
exercise  a  large  influence  on  the  flow  field. 

To  the  authors'  knowledge,  no  mathematical  solution  exists  which  will 
accurately  predict  the  flow  field  over  the  delta  wing  with  leading  edge  separation. 

Lift  Curve 

Figure  4  shows  the  I  ift  curves  of  the  delta  wing  models  tested  to  be  non- 

I  inear.   The  actual  lift  coefficients  exceed  those  predicted  by  slender  wing  theory, 

13,4 
It  has  been  theorized     that  the  extra  lift  produced  by  a  delta  wing  with  leading 

edge  separation  is  due  to  the  extra  entrainment  of  air  by  the  leading  edge  vortices, 

Furthermore,  that  the  non-1  inear  nature  of  this  extra  I  ift  is  due  to  the  inward 

and  upward  movement  of  the  vortex  cores  and  their  increase  in  strength  with  angle 

of  attack,  giving  an  ever  increasing  entrainment  effect. 

10  14 

Adams   and  Edwards   have  derived  formulas  for  predicting  the  lift  curve 

based  on  Legendre's  theory.  Their  formulas  are  applicable  only  for  small  angles 

of  attack.   Figure  4  shows  the  results  of  these  formulas  as  well  as  comparisons 


with  slender  wing  theory  and  experiments.  The  predicted  lift  curves  fall  short 
of  desired  accuracy. 
Adams'  formula: 


CL=  T[-fi!p<  +  77(1+   27r)(4A03  =* 


Edwards'  formula: 

cL=  f^°<  +  rr(\  +  27rH/^f  «y 

An  empirical  lift  curve  formula  has  been  developed  from  this  investigation 
which  is  covered  in  a  subsequent  section. 


MODELS  AND  APPARATUS 
Mode  I s 

All  models  were  aspect  ratio  one  delta  wings.  Model  I  was  a  flat  plate  of 
1/8  inch  thick  sheet  metal.  A  hollow  version  of  model  one  was  also  constructed 
for  the  smoke  flow  work.   See  Fig.  I. 

Models  II  and  III  were  Flberglas  models  of  thicknesses  12  percent  and  6  percent 
respectively.   The  root  chord  sections  are  parabolic  with  maximum  thickness  at 
the  0.50  root  chord  location,  as  shown  in  Fig.  I. 

For  smoke  flow  studies,  holes  of  1/16  inch  diameter  were  placed  at 
appropriate  locations  on  the  upper  surfaces. 

Smoke  Generator 

Kerosene  is  used  for  producing  the  smoke.  The  kerosene  is  boiled  by  passing 
it  through  an  electrically  heated  tube.  The  vapor  passes  through  a  small  orifice 
and  is  broken  into  tiny  particles,  which  produces  a  dry,  dense  white  smoke. 

W  ind  Tunnel s 

The  Princeton  4  ft.  by  5  ft.  wind  tunnel  was  used  for  force  measurements  and 
the  Princeton  3  ft.  by  4  ft.  three-dimensional  smoke  tunnel  was  used  in  the  smoke 
flow  visualization  studies. 
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PROCEDURE 
Force  Measurements 

The  I  iff  force  was  measured  by  the  beam  balance  system  of  the  wind  tunnel. 
The  tunnel  velocity  was  115  fps,  giving  a  Reynolds  Number  based  on  root  chord  of 
approximately  1.5  x  10  . 

Smoke  Flow  Visualization 

The  technique  can  best  be  described  by  a  three-dimensional  view  of  the  smoke 
tunnel  as  shown  in  Fig.  5.  The  model  is  mounted  in  the  test  section  as  shown. 
For  the  vortex  core  pictures  in  the  cross-flow  plane,  a  1600  mm  projector  with  a 
1000  watt  bulb  was  used  to  project  a  plane  of  I  ight  normal  to  the  free  stream  as 
shown  on  Fig.  6b.   The  plane  was  created  by  inserting  a  piece  of  metal,  having  a 
1/32  inch  slit  in  front  of  the  bulb  casing  as  shown  in  Fig.  5b.   Smoke,  under 
pressure,  is  forced  through  a  tube  into  the  hollow  model  and  out  the  holes  in  the 
wing.   Tunnel  velocity  was  maintained  at  II  fps.   The  photographs  of  the  vortex 
cores  were  taken  through  a  hole  in  the  downstream  section  as  shown.   By  means  of 
this  procedure,  quantitative  data  was  obtained  of  the  core  positions. 

The  distance  from  camera  to  model  was  approximately  twenty  feet.   The 
pictures  were  taken  with  a  35  mm  camera  with  a  360  mm  telephoto  lens.   With 
f  =  5.5,  pictures  of  the  flow  on  the  model  were  taken  at  one  second  exposures  and 
downstream  shots  were  taken  with  two  second  exposures. 

For  the  top  and  bottom  views  of  the  flat  plate  model,  a  Speed  Graf  I  ex 
camera  was  used.   Light  was  supplied  by  standard  flood  lights. 


9. 
RESULTS  AND  DISCUSSION 
Flow  Visual  ization  in  the  Cross-Flow  Plane 
Flat  Plate  Delta 

Figures  7  and  8  show  the  vortices  at  0.50  root  chord  and  the  trailing 
edge,  respectively,  at  various  angles  of  attack.   The  inward  and  upward  movement 
of  the  cores  as  angle  of  attack  is  increased  is  evident,, 

Figure  9  shows  the  vortices  at  0.25  cr  downstream  of  the  trailing  edge. 

o 
Up  to  ck  =  20,  the  primary  and  secondary  vortices  appear  to  be  comparable  in 

size.   Above  0<  =  20°,  the  secondary  vortices  become  smaller  than  the  primary 

vortices.   This  change  in  relative  size  may  indicate  a  reduction  in  the  effect 

of  the  secondary  vortices  on  the  flow* field  as  angle  of  attack  is  increased. 

Figure  10  shows  a  sequence  at  various  chordwise  stations  on  the  model  at 
C<   =  20  .   The  build  up  of  the  vortices  as  the  flow  moves  toward  the  trailing 
edge  is  shown. 

The  relation  between  the  primary  and  secondary  vortices  and  the  feeding 
sheet  at  locations  from  the  trail  ing  edge  to  0.50  cr  downstream  of  the  trailing 
edge  for  C<  =  15  ,  is  illustrated  in  Fig.  II.   Since  the  photographs  do  not 
show  clearly  the  position  of  the  feeding  sheet  at  the  various  downstream  locations, 
the  sketches  have  been  included.   At  the  trailing  edge,  the  triangular  regions 
outboard  and  below  the  primary  vortices  are  the  viscous  regions  in  which  the 
secondary  vortices  have  formed.   At  0.10  cr  downstream,  the  wake  intersects  the 
sheet  between  the  primary  and  the  secondary  cores.   This  point  of  intersection  is 
a  stagnation  point  in  the  cross-flow  plane  and  the  flow  sp I  its  and  proceeds  toward 
the  primary  and  the  secondary  cores.   Because  this  stagnation  point  occurs  on  the 
sheet  between  the  two  vortices,  it  indicates  that  the  effects  of  the  pressure 
gradients  produced  by  the  two  vortex  systems  are  almost  equal  at  that  point.   From 
this,  it  appears  that  the  strength  of  the  secondary  vortices  is  significant 
compared  to  that  of  the  primary  vortices.   Figure  27  shows  a  total  head  survey 
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made  in  the  cross-flow  plane  at  0.67  cr  on  a  model  similar  to  the  flat  plate  delta 
tested  in  this  report.   The  smoke  patterns  shown  in  Fig.  II  are  quite  similar  to 
the  contours  found  by  the  survey,   Figure  II  shows  that  as  the  flow  proceeds 
downstream,  the  feeding  sheet  moves  toward  the  primary  vortex,  and  the  secondary 
vortices  rotate  about  the  primary  vortices. 
6%   Delta 

Figure  12  shows  a  sequence  at  various  chordwise  stations  on  the  model 
at  0<  =  15  .   Figure  13  shows  the  vortices  at  the  trailing  edge  at  various  angles 
of  attack.   Qualitative  comparison  of  the  flat  plate  delta  and  the  6%   delta  is 

difficult.   However  it  does  appear,  from  a  comparison  of  Fig.  8a  and  Fig.  I3a,that 

o 
at  CX  =  5  ,  .  the  vortices  at  the  trailing  edge  of  the  flat  plate  delta  are  well 

defined  compared  to  the  vortices  on  the  6%   delta. 

Figure  14  shows  the  cross-flow  plane  at  0.25  cr  downstream  of  the  trailing 

edge.   The  position  of  the  feeding  sheet  between  the  primary  and  secondary  vortices 

is  shown  clearly  in  Figs.  14c  and  I4d.   The  diameters  of  the  cores  and  the 

vortices  increase  with  angle  of  attack.   It  appears  that  the  relative  size  of  the 

secondary  and  the  primary  vortex  is  not  as  significant  for  the  6%   delta  as  it  was 

for  the  flat  plate  delta.   For  example:  Fig.  9b  (flat  plate  delta,  OC=  10  ), 

shows  the  secondary  vortex  to  be  almost  equal  in  size  to  the  primary  vortex  and 

Fig.  14b  (6%   delta,  CX.  =  10  ),  shows  the  secondary  vortex  to  be  approximately  one 

half  the  size  of  the  primary  vortex.   It  may  be  speculated  that  the  effect  of  the 

secondary  vortices  on  the  flow  field  diminishes  with  wing  thickness. 

Flow  Visualization  on  the  Surfaces  of  the  Flat  Plate  Delta 

It  can  be  seen  in  Fig.  15  that  at  CX  =  0,  there  is  laminar  flow  over  the 
wing.   At  this  angle  of  attack  the  lift  curve  slope  may  be  computed  from 
linearized  potential  theory.  This  has  been  done  with  good  accuracy  by  many 
invest  igators. 
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As  CX  is  increased  to  2  ,  it  can  be  seen  in  Fig.  16  that  the  flow  is 
beginning  to  form  a  vortex  along  the  leading  edge  on  the  upper  surface.   The 
under  surface  flow  is  relatively  unaffected. 

Turning  to  Fig.  17  for  CX  =  4  ,  the  spiral  vortex  along  the  leading  edge 
can  now  be  seen  more  clearly.   Very  little  effect  is  noted  on  the  under  surface, 
as  the  tendency  to  flow  toward  the  leading  edge  is  slight. 

In  Fig.  18  for  CX  =  6°,  the  spiral  vortex  continues  to  grow;  and  as  more  air 
is  entrained  in  the  vortex  it  is  seen  that  the  region  of  reattached  flow  is 
decreasing.   On  the  under  surface,  more  air  is  flowing  toward  and  over  the  leading 
edge  and  is  entrained  in  the  vortex. 

As  CX  increases  to  10  ,  the  vortex  grows  and  the  reattached  region  shrinks. 
On  the  under  surface,  the  amount  of  flow  that  is  passing  over  the  leading  edge 
and  being  entrained  is  increasing  slowly. 

In  Figs.  20  to  24,  the  flow  on  the  upper  and  lower  surfaces  is  shown  with 
increasing  angles  of  attack,  up  to  35  .   It  can  be  seen  that  more  of  the  upper 
surface  flow  continues  to  roll  into  the  vortex  and  the  reattached  region 
decreases  in  size? until  at  CX  =  35  there  is  essentially  no  reattachment. 

On  the  under  surface,  more  and  more  of  the  flow  is  passing  out  and  over  the 
leading  edge.   This  effect  is  very  pronounced  at  the  higher  angles  of  attack. 
It  may  also  be  seen  in  Fig.  25a,  which  Is  a  view  of  the  flow  over  the  leading 
edge,  and  also  in  Fig.  25b,  which  is  a  view  from  the  trailing  edge. 

Figure  26  shows  the  vortex  cores  downstream  of  the  trailing  edge  as 
viewed  from  the  top  and  is  an  illustration  of  their  tenacity. 
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Graphical  Analysis  of  Vortex  Core  Loci  on  the  Body 

Flat  Plate  Delta 

Examining  Fig.  28  it  may  be  seen  that  the  curves  of  z/a  and  y/a  fall 
on  one  line  for  the  chordwise  positions  chosen.   This  chordwise  independence 
is  also  illustrated  in  Figs.  31  and  32.   The  core  is  seen  to  be  independent  of 
chordwise  location  until  about  0.80  cr  where  the  core  moves  inward  and  upward. 
This  can  be  attributed  to  trailing  edge  effects;  that  is,  the  core  line  is 
bent  toward  the  direction  of  the  free  stream  as  it  nears  the  trailing  edge. 

The  core  positions  with  respect  to  z/a  and  y/a  versus  ex  /  £      are  plotted 
in  Figs.  29  and  30.   Here  it  may  be  seen  that  z/a  is  approximately  I  inear  with 
angle  of  attack  for  CX  >  £   ,   |t  is  also  seen  that  the  core  moves  inward  very 
rapidly  at  small  angles  of  attack.   At  CX  =  6  ,  the  core  has  already  moved 
inboard  about  90  percent  of  its  total  inward  movement.   Theoretical  curves  of 

Q 

Brown  and  Michael   have  been  included  for  comparison.   Their  z/a  versus  o</£" 
curve  shows  good  agreement  with  the  experimental  results. 

Figures  33  and  34  show  the  side  view  and  plan  view  with  the  z  and  y  scales 
amplified  to  illustrate  the  core  movements. 

6%   Delta 

Now,  turning  to  Fig.  35,  it  may  be  seen  that  z/a  versus  y/a  continues 
to  move  inward  and  upward  with  increasing  angle  of  attack.   Individual  movements 
are  seen  in  Figs.  36  and  37.   The  movement  is  no  longer  independent  of  chordwise 
position.   The  core  locus  moves  inward  with  increasing  chordwise  position.   These 
effects  are  clearly  shown  in  Figs.  38  and  39.   Again,  amplified  scales  are  used 
to  give  the  side  and  plan  views  of  the  core  in  Figs.  40  and  41. 

\2%   Delta 

In  similar  fashion  to  the  flat  plate  and  the  6%   delta,  the  \2%   delta 
exhibits  inward  and  upward  movement  with  increasing  CX  as  seen  on  Figs.  42,  43, 
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and  44.   Since  the  thickness  is  a  function  of  chordwise  location,  the  core 
position  is  not  independent  of  x.   In  Fig.  45,  it  can  be  seen  that  z/a  decreases 
with  increasing  x  until  the  core  is  swept  up  by  the  trailing  edge  effects. 
Figure  46  shows  that  y/a  is  relatively  constant  with  x.   The  effect  of  body  shape 
can  also  be  seen  in  Fig.  47,  where  the  curvature  of  the  z  versus  x  curves  is 
distinctly  related  to  the  increase  and  decrease  of  body  thickness  with  x. 
Figure  48  is  a  plan  view  showing  y  versus  x,  with  an  amplified  y  scale. 
Downstream  Core  Loci 

Because  of  the  tenacity  of  the  vortex  cores  behind  the  delta  wing  (Fig.  26), 
it  was  possible  to  obtain  core  position  data  as  far  downstream  as  1.5  chord 
lengths.   The  cores  were  observed  very  clearly  as  far  downstream  as  eight  chord 
lengths  (16  ft.). 

Flat  Plate  Delta 

In  Fig.  49,  the  variation  of  z'/jr  with  downstream  station  is  shown.   It 
is  immediately  observable  that  at  all  angles  of  attack^ the  z'  position  varies  in  a 
sinusoidal  manner  as  the  core  proceeds  downstream.   As  the  angle  of  attack 
increases,  the  path  of  a  vortex  core  slopes  downward  at  an  increasing  angle.   It 
is  seen  that  with  all  the  core  paths  for  the  various  angles  of  attack  plotted  on 
one  figure,  the  paths  cross  between  0.25  and  0.50  cr  downstream. 

In  Fig.  52,  the  variation  of  y/:r  with  downstream  station  is  also  seen  to  be 
of  a  sinusoidal  nature.   There  is  a  trend  toward  inward  motion  of  the  cores  with 
angle  of  attack.   This  is  just  a  continuation  of  the  inward  motion  of  the  cores 
over  the  body  as  ex  increases. 

6%   Delta 

In  Fig.  50,  the  z' /—  versus  %   cr  downstream  curves  possess  the  same 
characteristics  as  the  corresponding  flat  plate  curves,  however,  the  sinusoidal 
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variation  is  less  marked.   It  can  be  seen  that  the  paths  of  one  core,  for 
increasing  angles  of  attack,  cross  between  0.50  and  0.75  cr  downstream. 

In  Fig.  53,  the  y/|  versus  %   cr  downstream  curves  show  a  very  pronounced 
sinusoidal  relation.   The  core  loci  also  show  an  inward  motion  with  angle  of 
attack. 

\2%   Delta 

_b 
In  Fig.  51,  the  z'/o  versus  %   cr  downstream  curves  again  show  the 

same  characteristics  as  the  corresponding  flat  plate  curves.   A  slight  sinusoidal 

,b 
variation  is  observed.   The  y/2  versus  %   cr  downstream  is  also  sinusoidal^  as 

seen  in  Fig.  54.   The  paths  of  one  core,  for  increasing  angles  of  attack,  cross 

between  0.50  and  0.90  cr  downstream. 

Summary  of  Downstream  Effects 

The  vortex  cores  are  seen  to  spiral  downstream  on  a  he  I  ical  path  as 

shown  in  Fig.  55.   This  accounts  for  the  apparent  sinusoidal  variation  of  the 

b      b 

z'/^  and  y/2  coordinates  versus  downstream  location.   This  spiral  is  of  the 

same  sense  as  the  spiral  flow  around  the  core.   The  spiral  nature  of  this  core 
path  is  related  to  the  presence  of  the  secondary  vortices.   The  secondary  vortex 
is  not  negligibly  small  compared  to  the  primary  vortex,  and  thus  the  center  of 
gravity  of  the  vortex  system  is  not  at  the  center  of  the  primary  corej as  seen 
in  Fig.  52b.   It  is  known  that  the  secondary  vortex  rolls  up  around  the  primary 
core,  thus  it  is  possible  to  explain  this  motion  by  assuming  both  vortices  rotate 
around  the  center  of  gravity.   This  would  explain  the  rolling  up  of  the  secondary 
vortex  around  the  primary  and  the  movement  of  the  primary  core  on  a  hel  ical  path 
downstream. 

It  is  also  seen  that  the  paths  of  one  core  pass  through  the  same  downstream 
vertical  position  as  the  wing  moves  through  the  angle  of  attack  range.   This 
point  appears  to  move  downstream  slightly  as  the  wing  thickness  increases. 


Thickness  Effects  on  Core  Positions  over  the  Body 

The  effect  of  thickness  on  the  vortex  core  position  may  be  seen  in  Figs.  56, 
57,  and  58, which  are  plots  of  core  positions  at  0.50  cr.   It  is  seen  that  as 
thickness  increases,  the  core  moves  outward  in  an  approximately  linear  manner. 
It  may  also  be  seen  that  increases  in  thickness  cause  no  noticeable  trend  of 
movement  of  the  core  in  the  vertical  direction. 
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EMPIRICAL  APPROACH  TO  THE  LIFT  CURVE  FOR  LOW  ASPECT  RATIO  DELTA  WINGS 


It  is  weil  known  that  the  C[_  versus  (y  curve  for  low  aspect  ratio  delta  wings 
has  the  shape  shown  in  Fig.  59.  That  is,  the  C[_  exceeds  that  predicted  by 
linearized  potential  theory  due  to  the  vortices  formed  by  leading  edge  separation. 

First,  let  us  assume  that  the  I  Iff  curve  is  composed  of  a  linear  component 
and  a  A  C[_  due  to  the  vortices  as  shown  in  Fig.  60. 


where 


and 


that  is: 


(I)   CLS  CL,+  AC, 


2TTA? 
C.  =  p/p  l;  CX.    (R.  T.  Jones'  expression  for  low  aspect  ratio 

*  delta  wings) 


A.  Ci  =  f  (vortex  position,  vortex  strength) 
CL  =  f  (y/a,  z/a,T  ) 


Expressing  in  total  differential  form  and  linearizing! 

m    dACu,|^dVQ  +  |^d%  +f^dr 

Curves  of  A  C\_   versus  z/a  and  y/a  have  been  plotted  in  Figs.  61  and  62.   To 
evaluate  the  partial  derivatives,  one  must  obtain 

(MCL\  pACj,  /BACA 

V  d  Y/a  )z/ot  p-conrt.    j  {BVci     l^/aT-  Const.  \    3T      /Y/q,2/a:Con5t 


ar_ 


To  do  this,  assume  ^—  =  constant.   This  approximation  is  based  on  the 
theoretical  calculations  of  V   versus  cX  in  Ref.  8.   Now,  any  non-1  inearity  in  the 
I   versus  oc  curve  is  a  second  order  effect  and  will  be  neglected. 

With  this  assumption,  it  is  possible  to  evaluate  the  first  two  derivatives  in 
equation  (2)  by  hoi  ding cx= constant  instead  of P  =  constant. 

Examining  Figs.  61  and  62,  it  can  be  seen  for  a  typical  case: 

« -18* :  V*  ■  030     (F,g  6i)  ;  &&  =  &2S  ,  0.  e 73  (Pf ,  e2) 
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cx-l8°:  Vo=  0.716  (f,-<5.62);  T^-^§f§  =  '-6&  (f^.ii) 


These  partial  derivatives  differ  by  a  factor  of  2.5,  so,  as  a  first  approximation 

it  will  be  assumed  that  x  7  .L  ^>/   g;  V/L      ,  and  thus  x-tt,l    will  be  neglected. 

a  va      o  T/o  $  '/a 

The  derivative  -^-=L  must  now  be  evaluated.   In  Ref.  15,  it  is  stated  that  the 
distance  between  the  vortex  cores,  "2Y",  is  a  measure  of  the  vortex  strength, 

(3)      2V  =  _pvr 

It    is   readily    seen   that   by   differentiating  equation    (3)    and   using 
,a>         aACu,   3ACL    dy/a 

that 

_^AC.  _   aACL  /-Y 

an       3%  voir, 

Since  it  has  been  assumed  that     _  ... L  <C^  .  7  .  L  un  less   ,_,  >^>  I,  the  term 

-.  p  u  will  be  neglected.   It  is  known  that  y/a  is  always  less  than  one  and  from 

Ref.  9  that  T»  I.   Therefore, 
^ACL  ,s    _1ACL 

and  will  be  neglected. 

In  view  of  the  above,  Eq.  (2)  may  be  reduced  to; 

where       *"  is  obtained  from  Fig.  61. 

The  equation  will  first  be  evaluated  for  the  flat  plate  case. 

The  relation  of  A  C|_  versus  z/a  was  obtained  by  Lagrange's  polynomial  formula. 

(5)  ACL  =  4-  (^o)Z- 0.  8  (5&)  +  0.0  6  . 

and      ^^  =-8%  -0,6  0 
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(6)  %  =  0.616  a  +  0.088 

Substituting,  the  equation  becomes; 

(7)  ACL=  1.52  cxz-0,0594  ex 

and  the  total  expression  becomes, for  A^  =  I,  t/c  =  0: 

(8)  CL=  1.48*  +-  L52o<2 

Having  an  expression  for  C|_  for  /£  =  I  and  t/c  =  0,  it  is  possible  to  expand  this  to 
other  thicknesses  and  aspect  ratios. 

For  aspect  ratios  other  than  I,  with  t/c  =  0: 

o)     ACL  =  0.0^25  fye)*-  0.0 14-6  <%= 

or 

(9a)    ACL=  0.0^25  (^t-J  -0.0146  ^rf 

Using  the  experimental  data  for  the  thickness  effect  it  can  be  derived  from  Fig.  63 
that: 

(io)      AC  ^  -  (0.  52^^-  0.034)  i% 

Combining  the  expressions,  an  equation  is  obtained  which  should  be  applicable 
over  a  range  of  smal I  /^  's  and  t/c' s. 

These  equations  were  then  appl  ied  to  delta  wings  of  aspect  ratios  of  one  to  two 
and  t/c  from  0  to  12  percent. 

The  results  are  shown  in  Figs.  64,  65,  66  and  67.   It  is  seen  that  the 
equation  gives  very  good  results.   From  this  it  can  be  concluded  that  the  additional 
lift  due  to  the  vortices  is  due  primarily  to  the  vertical  height  of  the  vortex  cores 
above  the  plane  of  the  wing. 

To  simplify  calculations,  equation  (I)  can  be  expressed  as: 

CL=A+B-C-D  where 
a  _  2-rr/^  or  *  <->         °"°  3 25  cX2"   '  r>  OC  /  \JSL 


A,  B,  C  and  D  are  plotted  on  Figs.  68,  69,  70  and  71  for  convenience  in 
evaluating  a  particular  case. 


20 


CONCLUSIONS 

Smoke  flow  visualization  techniques  can  be  a  powerful  tool  in  the 
investigation  of  flow  over  a  delta  wing. 

The  path  of  the  vortex  cores  downstream  of  the  trailing  edge  is  a  helix. 
This  helical  motion  is  caused  by  the  rotation  of  the  secondary  vortices  about 
the  primary  vortices. 

The  strength  of  the  secondary  vortices  relative  to  that  of  the  primary 
vortices  is  not  negligible.   Hence,  no  mathematical  treatment  which  fails  to 
include  the  effects  of  the  secondary  vortices  on  the  flow  field  can  accurately 
predict  the  I  i f t  curve.  , 

The  lift  may  be  considered  to  consist  of  a  linear  component,  wh ich  can  be 
predicted  by  potential  theory,  and  an  additional  lift  due  to  leading  edge  vortices, 
An  empirical  expression  has  been  developed  for  this  additional  lift,  which 
relates  it  to  the  height  of  the  primary  vortex  cores  above  the  wing.   This 
expression  predicts  the  lift  curves  accurately  for  sharp  leading  edged  delta 
wings  of  various  thicknesses  and  of  aspect  ratios  from  one  to  two. 
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